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We formulate equations of motion and conservation laws for a quantum many-body system in a 
co-moving Lagrangian reference frame. It is shown that generalized inertia forces in the co-moving 
frame are described by Green's deformation tensor g^vijsji) and a skew-symmetric vorticity tensor 
Ffi, v (£,t), where £ in the Lagrangian coordinate. Equations of motion are equivalent to those for a 
quantum many-body system in a space with time-dependent metric g^v(£,t) in the presence of an 
effective magnetic field F^, u (^,t). To illustrate the general formalism we apply it to the proof of the 
harmonic potential theorem. As another example of application we consider a fast long wavelength 
dynamics of a Fermi system in the dynamic Hartree approximation. In this case the kinetic equation 
in the Lagrangian frame can be solved explicitly. This allows us to formulate the description of a 
Fermi gas in terms of an effective nonlinear elasticity theory. We also discuss a relation of our results 
to time-dependent density functional theory. 
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I. INTRODUCTION 

Lagrangian and Eulcrian formulations of fluid mechan- 
ics are known as two alternative ways to describe dynam- 
ics of continuum media 1 . The more common Eulerian (or 
spatial) formulation considers basic collective variables, 
such as density n(x, t) or current j(x, t) distributions, as 
functions of space-time coordinates x and t ' . This cor- 
responds to the description of a system from the standard 
point of view of an observer in a fixed laboratory refer- 
ence frame. Central notions of Lagrangian (or material) 
description are the trajectories of infinitesimal fluid ele- 
ments. Every small element of a fluid can be uniquely 
labeled by its initial position £ that plays a role of inde- 
pendent, so called Lagrangian, coordinate. Lagrangian 
description represents the dynamics of continuum media 
as it is seen by a local observer, moving with a flow. In 
the last decades the Lagrangian method attracts an in- 
creasing attention as a powerful tool for studying nonlin- 
ear dynamics of compressible media with numerous ap- 
plications in cosmology, plasma physics, physics of semi- 
conductors, etc. (for a recent comprehensive review see 
Ref. 3). Recently we have shown that the Lagrangian 
coordinate naturally appears in time-dependent density 
functional theory (TDDFT), where it plays a role of a 
basic variable for a nonadiabatic exchange correlation 
potential 4 . It is also interesting to note a relation of 
Lagrangian fluid dynamics to noncommutative geometry 
and noncommuting gauge fields 5 . 

Commonly Lagrangian and Eulerian descriptions are 
considered as inherent ingredients of the classical con- 
tinuum mechanics. In fact, they offer two alternative 
techniques for solving the equations of classical hydrody- 
namics. However, the main idea of Lagrangian method, 
which is the description of dynamics using co-moving co- 
ordinates, is clearly much more general and universal. In 



the present paper we formulate microscopic equations of 
many-body dynamics in the co-moving Lagrangian refer- 
ence frame. The transformation to the Lagrangian frame 
corresponds to an explicit separation of the convective 
motion of particles. This is a natural generalization of 
the common separation of the center-of mass motion in 
homogeneous many-body systems. The separation of the 
center-of-mass motion also plays an important role in the 
theory of harmonically trapped systems. For the har- 
monic inhomogeneity the convective motion can be sep- 
arated by the transformation to a global accelerated ref- 
erence frame, which is a key step in the proof of the har- 
monic potential theorem 6,7,8 (HPT). In fact, the proof of 
HPT can be viewed as the simplest application of the 
Lagrangian description to quantum dynamics. In the 
case of a general inhomogeneous flow the separation of 
convective "center-of-mass" motion leads to an appear- 
ance of inhomogeneous inertia forces in the equations for 
the relative motion. We show that these forces can be 
uniquely described by the symmetric deformation tensor 
g^v{€,,t) and a skew-symmetric vorticity tensor F^(£, t). 
The deformation tensor enters equations of many-body 
dynamics as an effective time-dependent metric, while 
the vorticity tensor plays a role of an effective magnetic 
field. 

A great advantage of the Lagrangian description of 
many-body dynamics is that in the co-moving frame both 
the density of particles and the current density become 
the exact integrals of motion. The current density is zero 
in every point of Lagrangian £-space, while the parti- 
cles' density distribution preserves its initial form. These 
"conservation laws" are guaranteed by a fine local com- 
pensation of inertia forces, external forces, and the force 
of internal stresses. The above force balance follows the 
local momentum conservation law (the exact microscopic 
Navier-Stokes equation) after the transformation to the 
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Lagrangian frame. We explicitly demonstrate that the 
exact internal stress force takes a form of a covariant di- 
vergence of a symmetric second-rank stress tensor. As a 
byproduct of our formalism we obtain a microscopic rep- 
resentation for the local stress tensor in a general quan- 
tum many-body system. 

The concept of quantum stress has been introduced by 
Schrodingcr in 1927 9 . Over the last two decades there 
has been a growing interest in understanding properties 
of quantum systems, such as molecules or solids, in terms 
of the stress density (see, for example, Refs. 10,11,12,13 
and references therein) . A derivation of a microscopic ex- 
pression for the kinetic part of the stress tensor in quan- 
tum many-body system causes no problem. This sim- 
ple generalization of the one-particle result has been ob- 
tained in the classical paper by Martin and Schwinger 14 . 
However, the derivation of the microscopic form for the 
interaction related stress tensor turned out to be not that 
simple 11 ' 14,15 ' 16,17 ' 18 . In this paper we present two alter- 
native derivations of the symmetric form for the stress 
density, which has been obtained by Puff and Gillis in 
Rcf. 17. In particular, we show that this form is con- 
sistent with the definition of the stress tensor via the 
variational derivative of the energy with respect to the 
metric tensor. 

The structure of the paper is following. In Sec. II we 
consider the standard Eulerian form of the conservation 
laws in a quantum many-body system. In this section 
we also present a compact derivation of the microscopic 
expression for the exact stress tensor. Section III is de- 
voted to the formulation of quantum many-body theory 
in the co-moving Lagrangian frame. In Sec. IIIA the 
key notions of Lagrangian coordinate and of the defor- 
mation tensor are formally defined. The derivation of 
the equations of motion in an arbitrary local noninertial 
reference frame is presented in Sec. IIIB. Here we also 
derive the form of transformed many-body Hamiltonian 
and discuss the physical meaning of generalized inertia 
forces. In Sec. IIIC we derive local conservation laws, 
and present a complete formulation of the many-body 
problem in the Lagrangian frame. It is shown that this 
problem corresponds to the solution of the equations of 
motion for the relative motion, supplemented by the lo- 
cal force balance equation. The force balance equation 
plays a role of an additional gauge condition that fixes 
the reference frame. Section IV contains simple exam- 
ples of application of the general theory. In Sec. IVA 
we interpret the harmonic potential theorem 6 in terms 
of dynamics in the Lagrangian frame. In Sec. IVB we 
apply the general formalism to the study of semiclassical 
collisionless dynamics of a Fermi gas, and shortly discuss 
a connection of our approach to TDDFT. It is shown that 
in the regime of a fast long wavelength evolution the ki- 
netic equation in the Lagrangian frame can be solved 
explicitly In this case the behavior of the system is de- 
scribed by an effective nonlinear continuum mechanics, 
which, after the transformation to the laboratory frame, 
reduces to the generalized collisionless hydrodynamics of 



Refs. 19,20. In Sec. V we summarize our results. 



II. CONSERVATION LAWS IN THE 
LABORATORY REFERENCE FRAME: 
DEFINITION OF THE STRESS TENSOR 

In this paper we consider a system of N interacting 
particles in the presence of a time-dependent external po- 
tential [/ C xt (x, t) . The corresponding Hamiltonian takes 
the following standard form 



H=T+W+U 
f = 



(1) 

J ^t( x )^(x), (2) 
t?=i J dxd^io(|x-x'|)^t( x )^t( x ')^( x ')^( x ) } (3) 

U= /dx(7 oxt (x,t)V' t (x)^(x), (4) 



where w(x) is the interaction potential, and ift and if) 
are the field operators, which satisfy proper commutation 
relations 



[V t (x),V(x')]±=5(x-x'). 



(5) 



The upper (lower) sign in Eq. (5) corresponds to fermions 
(bosons), and [A, B]± = AB ± BA. Using Hamiltonian 
of Eqs. (l)-(4) we obtain Heiscnbcrg equations of motion 
for ^-operators 

d V 2 



+ 



J dx'w(\x-x'\)^{x')ip(x')ip(x) (6) 



Equation (6) allows to derive equations of motion for any 
physical observable as well as for any correlation function. 
The most important of these equations are the local con- 
servation laws or balance equations, which should be sat- 
isfied for an arbitrary evolution of the system. Below we 
concentrate on conservation laws for the number of par- 
ticles and for momentum. These local conservation laws 
follow the equations of motion for the density, n(x, t), 
and for the current, j(x, t) respectively. Computing the 
time derivative of the density operator, 



n(x,i) = V> f (x,t)V>(x,t), 



(7) 



we obtain the continuity equation that is the local bal- 
ance equation for the number of particles 



dn I 8 ^ = 

at dx» 

where n(x, t) = (n(x, t)) and 



(8) 
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Here the angle brackets stand for averaging with the ex- 
act many-body density matrix. Similarly using Eq. (6) 
we derive the equation of motion for the current, Eq. (9) 
(see, for example, Refs. 14,18) 

m ^ +F £n +F ^ +n _d_ Uext= ^ (10) 

Equation (10) has a meaning of the local force balance 
equation in the fixed laboratory reference frame. Vectors 
F^ ln and F™* in Eq. (10) correspond to the forces, which 
are related to the kinetic and the interaction effects re- 
spectively 



,kin_ 9 1 dip ^v 2 n)(ll) 



pint _ / 
" J 



dx v 2m \'dx^ dx v ' dx v dx^ 
,dw(\x — x'|) 



dx' 



dx» 



"P2(X,X ) 



(12) 



In Eq. (12) we introduced the notation p 2 (x, x') = 
(ijr (x)ip* (x')ip(x')ip(x)} for the two-particle density ma- 
trix. Obviously, the last term in the left hand side in 
Eq. (10) is the force produced by the external potential. 
The kinetic force of Eq. (11) has a form of a divergence of 
a symmetric second rank tensor. This automatically im- 
plies vanishing integral kinetic force, / F^ ln (x, t)dx = 0. 
The Newton's third law requires that the force F™ 1 of 
Eq. (12) should obey the same property, which is however 
by far not obvious. In fact, the possibility to represent 
Eq. (12) in a divergence form has been a subject of a long 
discussion in the literature 11 ' 14,16 ' 17,18 . An elegant sym- 
metric representation of the stress tensor has been pre- 
sented (unfortunately without derivation) by Puff and 
Gillis in Rcf. 17. Since this representation is of pri- 
mary importance for our paper, below we give a compact 
derivation of the Puff and Gillis result. 

The symmetry of the function p2{x, x') with respect 
to the permutation of coordinates allows us to transform 
vector Fj?*, Eq. (12), as follows 



(x)= J 
,dw(\x>\) 



,dw(\x — x'|) 



2 / ' /X dx" 1 



— [p 2 (x-x',x) +p 2 (x,x-x')] 



If ,dw(\x'\) , 

dx [p 2 (x + x,x) -p 2 (x,x-x)] 



dx 1. 



dw{\x'\) 
dx>» 



p 2 (x,x - x') 



where <9„ = Using an obvious operator identity 



«'v_ 1= f x ' V e Ax ' v 
Jo 



dX 



we arrive at the following final representation for the local 



force F}^ 



^" t(x) = 9^^ (X) 



(13) 



where W M „(x) is a stress tensor, which is responsible for 
the contribution of interparticle interaction to the force 
balance 21 



1 f x'^x'" 
W ^=-2j d *-W 



x' v dw{\x'\) 



d\x'\ 

x[ p 2 (x + \x',x- {1- \)x')d\ (14) 
Jo 

In the next section we will show that parameter A in 
Eq. (14) has a deep geometric meaning. It can be asso- 
ciated to the natural parameter for a geodesic (straight 
line in the present case) that connects two interacting 
particles (see also Appendix A). 

Equations (11), (13) and (14) show that the net in- 
ternal force, + i 7 ^ 11 *, is representablc in a form of 
divergence of a symmetric second-rank tensor H^. Ten- 
sor describes local internal stresses in the fluid. A 
contribution of the convective motion of particles to this 
tensor is known exactly 2 . It equals to the macroscopic 
momentum flow tensor, mnv^v^, where v = j/n is the 
fluid's velocity. It is convenient to separate this con- 
tribution explicitly and rewrite the conservation laws of 
Eqs. (8), (10) in the following familiar form 



d 

d d 



dx^ 



(15) 
(16) 



where D 



vV is the convective derivative and 



P^u = T^ v + Wp V is the exact stress tensor, which con- 
tains the kinetic, and the interaction, W^ v , contri- 
butions. The interaction stress tensor, W^, is given by 
Eq. (14), while the kinetic part, T^ u , is defined as follows 

T^ v = -^((qf,ip) 1 q u ip + {q v ip) ] %ip - ^^V 2 ^) (17) 

where q = — iV— mv is the operator of "relative" momen- 
tum which accounts for the above-mentioned separation 
of the macroscopic convective motion. 

Equations (15) and (16) form a basis for a hydrody- 
namic description of a nonequilibrim many-body sys- 
tem. According to the Runge-Gross mapping theorem 
of TDDFT 22 the exact many-body wave function/density 
matrix (for given initial conditions) is a unique functional 
of velocity v(x, t). Therefore the stress tensor P M „ is also 
a functional of v. Hence Eqs. (15) and (16) can be viewed 
as a formally closed system of equations that completely 
determine the dynamics of collective variables n(x, t) and 
v(x, t). This dynamics is governed by the external force, 
nd^Uext, and by the force of internal stress, d v P^ u . Since 
the convective motion has been explicitly separated from 
the stress tensor, only the relative motion of particles 
contributes to P^ v . A particular form of P^ v should be 
determined from the solution of a many-body problem 
in a reference frame moving with the "center-of-mass" 
velocity v(x, t). In the rest of the present paper we de- 
rive equations of many-body dynamics in this co-moving 
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frame and present simple illustrative examples of their 
solutions. 



III. QUANTUM DYNAMICS IN THE 
LAGRANGIAN FRAME 



A. Definition of the Lagrangian reference frame 

Co-moving or Lagrangian frame is a local noninertial 
reference frame which moves with the velocity v(x, t) of 
the fluid. Formally the transformation to the Lagrangian 
frame corresponds to a nonlinear change of variables x = 
x (£> which maps old coordinates x to new coordinates 
For a given velocity distribution the transformation 
function, x(£,i), is defined by the following initial value 
problem 



dt 



= v(x(£,t),t), x(£,0)=£ 



(18) 



If the function v(x, t) is continuous and satisfies the Lip- 
schitz condition in x, there exists a unique solution to the 
first order differential equation of Eq. (18) 23 . Therefore, 
under the above conditions on the velocity distribution, 
the map: x — > £ is unique and invertible. 

Physically the function x(£,t) corresponds to the tra- 
jectory of an infinitcsimally small fluid clement. Every 
fluid element (and therefore every trajectory) is uniquely 
labeled by the element's initial position - the Lagrangian 
coordinate £. The inverse function £ = £(x, t), which de- 
termines the transformation from the Lagrangian to the 
laboratory reference frame, recovers the initial position 
of a fluid element that at instant t arrives at the point x. 
The nonlinear transformation of coordinates, x = x(£, t), 



J 



induces a change of metric 
(dx) 2 = g^dC, 



dx a dx a 



(19) 



In classical continuum mechanics the symmetric second 
rank tensor g^ v (^,t), Eq. (19), is known as Green's de- 
formation tensor 1 . This tensor is normally used to char- 
acterize a deformed state of a system within the La- 
grangian description. The corresponding contravariant 
tensor, g^ v ', is defined as follows 



g^g*,, 



<5 M 



dx^dx" 



(20) 



Since the deformation tensor g^ has a meaning of the 
metric tensor in the Lagrangian £-space, it should play 
a key role in the description of many-body dynamics. It 
is quite natural to expect that the general equations of 
motion in the Lagrangian frame should reduce to those 
in a space with time-dependent metric 3^(£,i). Below 
we confirm this intuitive expectations by explicit calcu- 
lations. 

B. Equations of motion in a local noninertial 
reference frame 

In this subsection we derive quantum equations of mo- 
tion in a general noninertial (not necessarily Lagrangian) 
reference frame. The frame is defined by its velocity 
v(x, t), which enters the trajectory equation of Eq. (18), 
and thus provides a unique and invertible map: x — > 
As a first step in the derivation we perform a nonlinear 
transformation of coordinates, x = x(£,t), in the equa- 
tion of motion, Eq. (6), and in the commutation relations 
of Eq. (5) . Straightforward calculations lead to the result 



4 t m = 



1 1 d _ „„ d x d rr . 



m+ J dt'w(!te,)tf(z')WMt) (2i) 



r 



where U cxt (£,t) = U ex t(x(£,t),t), and field operators 
t) satisfy the following equal-time commutation re- 
lations 



[^(^(0]± = ^(£-0 



(22) 



To shorten the notations in Eq. (21) we omitted the ex- 
plicit time dependence in the argument of ^-operators. 
The first term in the brackets in Eq. (21) is the Laplace 
operator in a reference frame with metrics g^ v (see, for 
example, Ref. 24) , while the second term comes from the 
transformation of the time derivative in Eq. (6). This 
term is proportional to vector u^(£, t) that is the vector 



of velocity, transformed to a new frame: 

^«,t) = ||x(x(€, *),*)■ 



(23) 



The interparticle distance, in the argument of 

the interaction potential in Eq. (21) equals to a length 
of geodesic that connects points £ and Geodesic, 
z^^'(A), parameterized by a parameter A (0 < A < 1), is 
a solution to the following equation 24 



^(A)+r^(z)i Q (A)i' 3 (A)=0, 
where z = dz/dX, and is affine connection 

dg v p dg a f}\ 



1 a/3 



— Q 



(24) 



(25) 
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Equation (24) should be solved with boundary conditions 
z (0) = £> z (l) = It is convenient to parameterize 
geodesies by a natural parameter, which is chosen in such 
a way that an absolute value of the "velocity", |z| = 
^Jg^ v z^z v , becomes independent of A along the curve 
z(A). For this parameterization the length which 
enters Eq. (21), simply equals to |z| at any point on the 
geodesic 



(26) 



Equation (21) is the equation of motion for the op- 
erator V(£>*) = ^( x (£> *)>*)■ D uc to the Jacobian fac- 
tor in the commutation relations of Eq. (22), the 
quantity ip(£,t) can not be interpreted as an operator for 
annihilation of a particle in a given point of £-space. In 
particular, the operator n(£) = ^ (£)?/>(£) does not cor- 
respond to the density operator in the new frame. By 
definition the density is a number of particles per unit 
volume that is changed under a volume non-preserving 
coordinate transformation. Therefore it is natural to de- 
fine the physical field operators and the density operator 
as follows 

$(£)=9h(fi, ^(O=0M(€). (27) 

n(€) = ^ t (€W€) = Vff^ + (€W€). (28) 

which automatically accounts for the proper change of a 
unit volume in the deformed reference frame. Obviously 
the redefined field operators satisfy the common 

commutations relations 



(29) 



The renormalization of ^-operators, Eq. (27), is equiv- 
alent to the corresponding multiplicative redefinition of 
the many-body wave function. This redefinition is aimed 
to preserve the common probabilistic interpretation and 
the standard form of the normalization conditions in the 
new reference frame (similar arguments were suggested 
by Podolsky 25 in early days of quantum mechanics). 

Let us show that the renormalization of field opera- 
tors, Eq. (27), also simplifies the form of the equations 
of motion. First we note that the differential operator 
in the right hand side in Eq. (21) (first two terms in the 
square brackets) can be rearranged as follows 

_L 1 d_ ^_d_ d 



1 K^y/gK* . 

- m-Z- i 



y/9 2m 



v^V .1/9 



2^9 W 



yfW (30) 



where we introduced an operator of "kinematic" momen- 
tum in the noninertial reference frame: 



w ■ d 



(31) 



(Raising and lowering of tensor indices are performed 
according to the standard rules, i.e. — g liv v v or 
=g^K u .) 

Using the equation of trajectory x(£,i), Eq. (18), and 
the definition of metric tensor g^ u , Eq. (19), one can 
prove the following identity 



dg 4 1 d In g 
~~dT = 4 dt 



(- 



d 



(32) 



The quantity in the right hand side in Eq. (32) coincides 
with the last term in the right hand side in Eq. (30). 
Hence the sum of the corresponding term in the equa- 
tion of motion, Eq. (21), and of the time derivative of ip 
reduces to the following compact form 



d ^ , 1 ( 9 r~»\ 



ldgtip 
dt 



Substituting Eq. (30) into Eq. (21) and using Eq. (33), we 
obtain the final equation of motion for the renormalized 
field operator ip(£,t) 



' dt 



9 4 ^ 9 

y 2m y 



4 + Uext - m- 



+ 



(34) 



Equation (34) allows us to recover a form of the trans- 
formed Hamiltonian , tp], which, together with the 
commutation relations of Eq. (29), determines the dy- 
namics of the system 



H = T + W + U, (35) 
T = f dtsfg (K^g-^y |^ (K v g-*$) , (36) 

w=l[ Kdz'wiiteW^MHeMeMz), (37) 



U= I d£ ( f/ext - m-^- ) i>H 



(38) 



Equations (34)-(38) represent the main results of this 
subsection. Equation (34) is the Heisenberg equation of 
motion for the physical field operator, while Eqs. (35)- 
(38) establish the rules for the transformation of the 
many-body Hamiltonian to an arbitrary local noniner- 
tial reference frame. 

Formally the Hamiltonian of Eqs. (35)-(38) describes 
a system of quantum particles in the presence of an ef- 
fective vector potential mv(£, t) and an additional ef- 
fective scalar potential mv 2 /2. The particles live in a 
space with the time-dependent metric g^ui^,, t) and inter- 
act via pairwise potential which depends on the length of 
geodesic connecting a pair of particles. Additional "po- 
tentials" and a nontrivial metric tensor are responsible 
for generalized inertia forces exerted on a particle in a 
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general noninertial reference frame. To get a transpar- 
ent physical understanding of these forces it is instruc- 
tive to look on dynamics in the semiclassical approxi- 
mation. Since the most important inertial contributions 
enter only quadratic parts of the Hamiltonian (Eqs. (36) 
and (38)), we neglect for a moment the interaction, and 
consider an equation of motion for the Wigner function 



/p&*)= / - 



in a gas of noninteracting particles. In the semiclassical 
limit the Wigner function satisfies the following kinetic 
equation 



df P | dH(p,£) df p 
dt dp d£ 



dg(p,0 dfp 
d£ dp 



0, 



(39) 



where H(p, £) is the semiclassical Hamiltonian function, 
which corresponds the noninteracting part of Eq. (35): 



H{p,i)^'—{p l _ l - mv^){p v - mv v ) 



VnV 



+ L/ext - ni—< 



(40) 



Substituting Eq. (40) into Eq. (39) we get the result 



dt 



( dg af} PaP0 dv a 
V dip 2m ~ 'dip" 1 



dU e: 



dpu 



0. (41) 



Inertia forces do not explicitly show up in Eq. (41). The 
reason is that Eq. (41) is the equation for the function 
/ p which depends on the canonical momentum p. The 
physical velocity of a particle in the new reference frame 
is proportional to the kinematic momentum K = p — mv 
(i.e. dH/dp^ = K^/m). Therefore it is more natural 
physically to consider K as an independent variable in 
the kinetic equation. The distribution function of the 
kinematic momentum, /k(£ ; t), can be introduced as fol- 
lows 



/K(€,*) = / K+m v(€'*)- 



(42) 



Performing the corresponding change of variables in 
Eq. (41) we obtain the final semiclassical equation of 
motion for the distribution function f^($,t) in the lo- 
cal noninertial reference frame 



dt 



+ 



m d^" 
d 



dg a[i K a KP 
d£ v 2m 



[Text - m- 



dK v 



= 0, 



(43) 



where a skew-symmetric second rank tensor is de- 
fined as follows 



de 



dv v 

w 



(44) 



Since tensor F^ v vanishes for an irrotational flow, we 
name it the vorticity tensor 26 . In the next section 
Eq. (43) will be applied to the derivation of generalized 
collisionless hydrodynamics 19,20 . 

The expression in the square brackets in Eq. (43) con- 
tains all inertia forces. These are all the terms except for 
the external force, d^U ext - The first term in the square 
brackets is the linear acceleration force, while the last 
term is related to the kinetic energy of a moving frame. 
In a particular case of a homogeneously rotating frame 
the last term is responsible for the centrifugal force. The 
second and the third terms in the square brackets corre- 
spond to inertia forces that depend on a velocity of a par- 
ticular particle. The second term is the classical Coriolis 
force. This force is proportional to the skew-symmetric 
vorticity tensor, which defines a local angular velocity of 
the reference frame. The third, bilinear in particle's mo- 
mentum term is less common. The corresponding inertia 
force makes a free particle to move along a geodesic in 
a local noninertial frame. Indeed, the third term in the 
square brackets in Eq. (43) can be rewritten as follows 



2m d£" m af} 



(45) 



where we have used Eq. (25), which relates the affine 
connection to the metric tensor g^ u . The right hand 
side in Eq. (45) is easily recognized as a covariant com- 
ponent of the force in the equation of geodesic (see, for 
example, Eq. (24)). 



C. Conserving quantities and balance equations 

1. The continuity equation 

The first problem we address in this subsection is a 

proper definition of the current operator, jP, in a general 
noninertial reference frame. The easiest way to establish 

a form of is to derive the equation of motion for the 
density operator n(£, t) = t). Using Eq. (34) 

to compute the time derivative of the density operator 
we find that the desired equation indeed reduces to the 
common form of the continuity equation, 



dn dj^ 
1 — — = 

dt + d& ' 

if we define the current operator, j' i (^, t), as follows 



(46) 



(47) 



The standard form of the continuity equation, Eq. (46), 
should be considered as one more justification for the 
redefinition of field operators, Eq. (27). We would like 
to outline a very natural form of the current operator, 
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Eq. (47). Despite the presence of the Jacobian factors 
(y/g or g 1 / 4 ) in the Hamiltonian, they completely vanish 
in Eq. (47) (as well as in the definition of the density 
operator of Eq. (28)). 

From this point we restrict ourselves to the Lagrangian 
frame, which is the local reference frame, moving with the 
velocity v of the fluid. In this special case the continuity 
equation admits a very simple solution. Let us calculate 

the expectation value of the current operator j>, Eq. (47). 
This can be done, for example, by transforming the right 
hand side in Eq. (47) back to the laboratory frame, and 
by using Eq. (9) together with the definition of the ve- 
locity, v = j/n. The result takes an extremely simple 
form 



(48) 



Thus the current density is exactly zero in every point 
of the Lagrangian £-space. This is of course not surpris- 
ing, since an observer in the co-moving frame should not 
see any current. Combining Eq. (48) and the continuity 
equation of Eq. (46) we arrive at the conclusion that the 
density n(£, t) is independent of time 



»(£,*) = n(€,0) = no(£), 



(49) 



where no(x) is the initial density distribution. Therefore 
in the Lagrangian frame not only the number of particles 
N is an integral of motion, but the density itself is also 
a conserving quantity. Evolution of the density in the 
laboratory frame can be calculated with the following 
formula (see Eq. (28)) 



n(x, t) = 



n(£(x, t),t) 



n (£(x,t)) 



Vff(€(x,t),t) y/g(£(x:,t),t) 



(50) 



Equation (50) is, in fact, the explicit solution to the con- 
tinuity equation of Eq. (8), which defines the density 
n(x, t) as a functional of velocity v(x, t). 

Equations (48) and (49) demonstrate the main advan- 
tage of the Lagrangian frame for the description of many- 
body dynamics. In this very special reference frame the 
inertia forces are adjusted to get exactly zero current den- 
sity and therefore to keep the density of particles fixed 
during the whole evolution of the system. Equation (48) 
can be used to construct a complete many-body theory 
in the co-moving frame. The frame's velocity v enters 
the equation of motion, Eq. (34), as an external param- 
eter. Imposing the local "gauge" condition of Eq. (48) 
we specify the reference frame and thus get the complete 
theory with all quantities defined by the initial condi- 
tions. 



2. Local force balance in the Lagrangian frame 

Let us turn to the local momentum conservation law. 
In the laboratory reference frame it is given by Eq. (16) 
(or, equivalently, by Eq. (10)). Since in the Lagrangian 



frame the current density is zero, the local momentum 
conservation law should reduce to the zero force condi- 
tion - the inertia forces should exactly compensate the 
external force and the force of internal stresses. Below 
we derive an explicit form of this balance equation by 
the direct transformation of Eq. (16) to the Lagrangian 
coordinates £. First we express the vector of velocity v 
and the stress tensor P^ v in terms of the corresponding 
quantities, v and P^, in the Lagrangian frame 



dp 



vp{£,t), 



dx» 

_dpdx» p 
dx» dp ^ 1 



(51) 
(52) 



Equation (51) follows the definition of Eq. (23), while 
in Eq. (52) we adopted the standard tensor transfor- 
mation rules 24 . Substituting Eqs. (51) and (52) into 
Eq. (16), transforming the derivatives, and multiplying 
the result by dx^/dp, we obtain the following equation 



dx* d (dp _ 

mn dp at [d^ V0 



dx» dp d 
~dp dx v dp 
<9t/ ox t 



dp dx v ~ 7 
dx» dp $ 



dp 



= 



(53) 



The first term in Eq. (53) can be further simplified as 
follows 



dx» d dp, 
'dpdtdx^ 1 



vp 



dt 



vp 



dp dv^ 
'dx^'dp 



dt 



1 dvpv 13 
2 



dp 
(54) 

In the sequence of transformations in Eq. (54) we used 
an obvious identity §|s-§§f = £f , the trajectory equa- 
tion of Eq. (18), and the definition of v*, Eq. (23). A 
similar simplification of the second term in Eq. (53) is 
even more straightforward. One only needs to apply the 
chain rule for the calculation of the spatial derivative and 
take into account the following explicit representation for 
affine connection (see, for example, Ref. 24), 



dp d 2 x 



Q/3 dx-r dpdp ' 



(55) 



As a result we get a very natural expression for the second 
term in Eq. (53) 



dx^ dp d ( dp 3 dx v ~ 7 



dp dx v dp \dx» di 







p0 

a;p 



(56) 



where the semicolon is used to denote the covariant 
derivative. The covariant divergence of the stress ten- 
sor in Eq. (56) is defined as follows 24,27 



dP v 

pv _ P_ I yi> pa 

1 dy/gP" 



pa 
1 Lia r v 



V9 dp 



1 dffq/3 pap 

2 dp 



(57) 
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Substitution of Eqs. (54) and (56) into Eq. (53) leads 
to the final form of the force balance equation in the 
Lagrangian frame 



dvf, d 



Uc 



+ ^P; ; , = 0. (58) 



A direct comparison of the force term in the kinetic equa- 
tion of Eq. (43) and the term in the square brackets in 
Eq. (58) shows that the later is exactly the sum of the 
external force and two inertia forces that are indepen- 
dent of particle's momentum. These three forces are bal- 
anced by the force of internal stresses (the second term 
in Eq. (58)). The net force, exerted on every fluid ele- 
ment in the Lagrangian space, is zero, which results in 
a zero current density and a stationary particles' density 
distribution. It should be noted that the rest of inertia 
forces (those, which arc different for different particles in 
a fluid element) implicitly present in the kinetic part of 
the stress tensor P^. To see this more clearly we need 
to derive an explicit microscopic representation for this 
tensor. 



3. Stress tensor in the Lagrangian frame 

In general both kinetic, and interaction, W^ u , con- 
tributions to the total stress tensor, P^ v = T^ v + W M1/ , 
can be found using the common transformations rules 24 

T^,t) = ——T a ^,t),t), (59) 



(60) 



Here stress tensors, T a @(x,t) and W a( 3(x, t), in the labo- 
ratory reference frame are given by Eqs. (17) and (14) re- 
spectively. We shall however follow another route, which 
takes full advantage of geometric ideas we develop in 
this paper. The transformed many-body Hamiltonian 



of Eqs. (35)-(38), is an explicit functional of the met- 
ric tensor g^ v . Therefore we can find the required stress 
tensor by computing the variational derivative of the en- 
ergy with respect to the metric tensor. More precisely, 
we make use of the fact that under a small variation of 
the metric, the variations of the kinetic energy, Eq. (36), 
and of the energy of interparticle interaction, Eq. (37), 
are related to tensors T M „ and respectively 



<HT)= I d^5g>"'T„„ 



d^Sg^ (61) 



<*«"-> = / d^Sg^W^ = - f d ^Sg^W^(62) 



Such a definition of the stress tensors is closely related 
to the common definition of the energy-momentum ten- 
sor in general relativity (see, for example, Ref. 27). The 
main advantage of this definition is that it automatically 
gives a symmetric form of the stress/energy- momentum 
tensors. Recently a very similar approach has been used 
to derive a microscopic expression for the stress tensor 
in the equilibrium quantum many-body system within 
the local density approximation 13 . The paper of Ref. 13 
also contains a general discussion of the above geometric 
definition of the stress tensors in the context of nonrcla- 
tivistic quantum mechanics. 

The variation of the Hamiltonian should be taken at 
constant ^-variables (since they satisfy the equations of 
motion) 13,27 and at constant velocity v. The kinetic en- 
ergy operator T of (36) contains the metric tensor only 
in a form of g^ v , ^[g and g~i . Noting that 



5 9 * = \g l 9^Sg^, 



we can easily compute the variation of T and then iden- 
tify the kinetic stress tensor using Eq. (61). The result 
of the calculations takes the form 



I 



If we evaluate the right hand side in Eq. (63) for Euclid- 
ian metric, g M „ = 5^, we immediately recover Eq. (17). 
Therefore the commonly used symmetric form of the ki- 
netic stress tensor T M „, Eq. (17), is in exact correspon- 
dence with the geometric definition of Eq. (61). 

Calculation of the variation SW, Eq. (62), is a little 
bit more involved. Since the interaction Hamiltonian of 
Eq. (37) depends on g^ v only via the length of geodesic, 



we have 

S(W) = \J dv d V '6l v ^^^p 2 (v,v') (64) 

where p 2 (^,'^ / ) = (V' t ( T ?)V'^( I 7')V'( T 7')V'( r 7))- The next 
step is to compute the variation of the functional 
ln,n' [sw] i Eq. (26). Let A in Eq. (26) be a natural pa- 
rameter for a geodesic in the space with "unperturbed" 
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metric g^ u (not the "full" metric g^ + Sg^). For this 
parameterization the variation of l v ri i [g^] takes the form 



51 



21 



■n,v' Jo 



d\6g^(z(\))z»(\)z»(\) 



21 



where z(A) = z T?jT? /(A) is the geodesic which connects 
points r\ and r\' . Substituting Eq. (65) into Eq. (64), 
and using the definition of Eq. (62) we get the following 
representation for the interaction part of the stress tensor 



lr),rj' dlrjrf' 



(66) 



clidcan metric 



9»v 



In this case I 



77,77' 



Let us evaluate the right hand side in Eq. (66) at Eu- £>(0) = 0). The initial value problem with the external 

\r] — rf'\ potential of Eq. (67) is exactly solvable, which is known 
as the harmonic potential theorem (HPT) 6 . It is also 
known that HPT is related to the covariance of the time- 
dependent Schodinger equation under the transformation 
to a global accelerated reference frame 7,8 . Therefore our 
formulation of the many-body problem should be per- 
fectly suited to the demonstration of HPT. 

Within the present Lagrangian formulation one needs 
to find a selfconsistent solution to the many-body equa- 
tion of motion, Eq. (34), and to the force balance equa- 
tion, Eq. (58). Let us assume that velocity v(x, i), which 
defines (via Eq. (18)) the motion of the reference frame, 
is a function of t only: v(x, t) = V(i). In this case the 
trajectory of a fluid element takes a form 



while the geodesic (parameterized by the natural param- 
eter) is a straight line 

z v,v'( x ) =V + (v' - »7)A 

The above expressions for and z ?7i77 '(A) should be 
substituted into Eq. (66). Introducing a new variable 

= jj' — rj, and removing the delta-function by the inte- 
gration over rj we obtain at the result that exactly coin- 
cides with Eq. (14). Therefore the symmetric representa- 
tion of Eq. (14), which has been obtained in the previous 
section by somewhat artificial manipulations, has a clear 
geometric meaning. In particular, the internal parame- 
ter A in Eq. (14) is the natural parameter for a geodesic 
connecting two interacting particles. 

Equations (63) and (66) are the principal results of 
the present subsection. They define explicit microscopic 
representations for the stress tensors in a local noninertial 
reference frame. 

The zero force condition of Eq. (58) with of 
Eqs. (63) and (66) is equivalent to the requirement of 
zero current density, Eq. (48). Hence we can use Eq. (58) 
as an alternative "gauge" condition to fix the velocity pa- 
rameter v(£, t), entering many-body equations of motion, 
Eq. (34). 



IV. EXAMPLES AND APPLICATIONS 

A. The harmonic potential theorem 

As a first simple example of application of our gen- 
eral formalism, we consider a many-body dynamics in 
the presence of the following external potential 

E/ext(x, t) = ^mu^x" + E^{t)x» (67) 

is a constant tensor and E^it) is a time- 



x(£,t)=£ + R(t), (68) 
where R(i) is a solution to the following Cauchy problem 
dR(t) 



dt 



= V(i), R(0) = 



where 

dependent vector (without loss of generality we can put 



j Liu 



Clearly, if our anzatz, v(x, t) — V(i), is a selfconsistent 
solution, then R(t) should correspond to the center-of- 
mass coordinate. Using Eq. (68) we get the following 
results for the metric tensor g^ v , the velocity v(£, i), and 
the effective potential, which enter Eqs. (34) and (58) 

ff„„ = <W. v(£,t)=V(t) (69) 

Uext(x(€, *),*)- = ^muv^r - L{t) 

+e[moJ^R v (t)+E^t)} (70) 

Function L(t) in Eq. (70) is the classical Lagrangian for 
a particle moving in the harmonic potential of Eq. (67) 

L{t) = y V 2 (t) - ±mcj^W(t)R v (t) - E^t)R»(t) 

The equation of motion, Eq. (34), and the force balance 
equation, Eq. (58), simplify respectively as follows 
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2m 2 p 



= 



^' (71) 
(72) 



where we also performed a gauge transformation, 
= ^(£> t) ex P[ — imY^ + i J Q Ldt'], which corre- 
sponds to the transformation from the canonical to the 
kinematic momentum. 

Let initially the system is prepared in a stationary state 
(or in arbitrary mixture of stationary states). This means 
that at t = the stationary force balance equation is 
fulfilled 



() 



0. 



(73) 



If at all f > the center-of-mass coordinate R(t) satisfies 
the classical equation of motion 



d 2 R»{t) 
1 dt 2 



+ moj^R v (t) + E^(t) = 0, 



(74) 



then both the equation of motion, Eq. (71), and the force 
balance equation, Eq. (72), preserve their initial (sta- 
tionary) form. Therefore the many-body system in the 
co-moving frame remains in the initial stationary state, 
in particular, P^ v (^,t) = P M „(£,0). This statement is 
the essence of HPT 6 . Within the present formulation of 
many-body dynamics it appears quite naturally. In fact, 
HPT is a built-in property of our sclfconsistcnt approach. 
This actually means that any approximate treatment of 
a selfconsistent system of Eqs. (34), (58) should automat- 
ically satisfy HPT. 



B. Geometric formulation of generalized 
hydrodynamics: nonlinear elasticity of a collisionless 
Fermi gas 

The HPT type of motion provides an extremely simple 
example of many-body dynamics without any deforma- 
tion of local fluid elements {g^ T = 5^ u ). In this sub- 
section we apply our approach to a much more general 
situation with a nontrivial dynamics of a fluid. Namely, 
we consider a semiclassical dynamics of an interacting 
Fermi system in the time-dependent Hartree approxima- 
tion. The problem reduces to the selfconsistent solu- 
tion of a semiclassical collisionless kinetic equation (see 
Eq. (43)) 



dt 



rn 



dg a p K<*K? 
di v 2m 



+ 



0f 



U — m- 



dK v 



(75) 



and a force balance equation (see Eq. (58)) 



m- 



d_Ufj_ 

dt 



d 



U 



n 



0. (76) 



Here U — U eict (£,t) + Un(£,t) is a sum of the external 
potential and the Hartree potential, 

Efa(€,*) = / «>(*€,*' )no(€')de'- 

Since the interaction effects are already included (on the 
mean field level) in the selfconsistent potential, only ki- 
netic part of the stress tensor contributes to Eq. (76), i.e. 
Pixv = m^ 1 X)k K/iKv flc/ \/9- The last expression for 
P^ v is a plain semiclassical limit of the general kinetic 
stress tensor, Eq. (63). 

The problem of solving Eqs. (75), (76) can be refor- 
mulated as follows. Let us substitute the sum of the in- 
ertia and the external forces from the balance equation, 
Eq. (76), into the kinetic equation of Eq. (75). After the 
substitution the potential U in Eq. (75) cancels out and 
the kinetic equation reduces to the following universal 
form 



dt 



m d^ v 
dg af3 K a K@ 
d£ v 2m 



K^F, 



[IV 



n ^ 



dK v 



(77) 



where the stress tensor P^ u is defined as follows 

a.t)4E^K.t) (78) 
v9 K 

Equations (77) and (78) constitute a closed set, which 
is structurally similar to the system of Vlasov equations 
with a selfconsistent force. The skew-symmetric vorticity 
tensor, F^ v (^,t), and the symmetric deformation tensor, 
9p,v{$,,t)i enter Eqs. (77) and (78) as external parame- 
ters, which govern the evolution of the system. Hence 
these equations define a distribution function, 
as a unique functional of F^ and t^„, provided the ini- 
tial condition, /k(£>0) = /k^(£)i i s gi ven - Equation (78) 
determines the stress tensor as a universal (i. e. indepen- 
dent of external potential) functional of F^ and g^ v : 



P\1V \P\IV1 9flis] (£) t) ■ 



(79) 



11 



The vorticity and the deformation tensors contain nine 
independent scalar functions (three from and six 
from g^) which completely describe a deformed state 
of a system 1 . Hence Eq. (79) plays a role of a general- 
ized "equation of state" which relates the stress tensor 
to the deformation. It is worth mentioning that the exis- 
tence of such an equation of state is a direct consequence 
of Rungc-Gross mapping theorem 22 in TDDFT. 

Substituting the functional of Eq. (79) into Eq. (76) we 
obtain a hydrodynamic equation of motion which deter- 
mines the evolution of velocity for a given external poten- 
tial. Therefore the description of many-body dynamics 
consists of two separate problems. The first one corre- 
sponds to the universal kinetic problem of Eqs. (77), (78). 
By solving these equations we find the stress tensor func- 
tional, Eq. (79) (the generalized equation of state). The 
second problem is to compute the velocity and density 
distributions by solving the closed set of hydrodynamics 
equations, Eqs. (18) and (76). 

The universal kinetic problem of Eqs. (77), (78) can 
be solved explicitly in the case of a fast long wavelength 
dynamics, i. c. if the deformation tensor is a fast func- 
tion of time, but slowly changes in space. More precisely, 
we assume that the characteristic length scale, L, of the 
deformation inhomogeneity is much large than tu, where 
r is the time scale of a dynamical process and u is the 
characteristic velocity of a particle. This situation is, for 
example, common in Coulomb systems where the plasma 
frequency determines the characteristic time scale of dy- 
namics, while the corresponding spatial variations of the 
density can be arbitrary slow. Let us estimate different 
terms in Eq. (77) under the above assumption. The first 
and the second terms in the right hand side in Eq. (77) 
are of the order of 1/r and u/L respectively. Both terms 
in the second line in Eq. (77) also give a contribution 
~ u/L, while the term, related to the Coriolis force, is 
proportional to F ~ Vt/L, where vt is a rotational (or 
transverse) component of the velocity. According to the 
force balance equation of Eq. (76), for any physical veloc- 
ity the transverse part of the linear acceleration is com- 
pensated by the transverse part of the vector ^P^. v ■ 
Hence vt should be proportional to tu 2 /L. This means 
that the contribution of Coriolis force to the kinetic equa- 
tion is of the order of tu 2 /L 2 . Therefore, to the leading 
order in the small parameter 7 = tu/L <C 1, only the 
first term in Eq. (77) gives a nonvanishing contribution. 
Thus the universal problem of Eqs. (77), (78) reduces to 
the following trivial equation 



|/k(^)=0. 



(80) 



Equation (80) shows that for a fast, small-gradient evo- 
lution the distribution function in the Lagrangian frame 

preserves its initial form: fj^(^,t) = /k^(£)- ^ n this 
respect the dynamics reminds the HPT type of motion. 
However the evolution of the velocity is by far not trivial. 
Below we consider a system which evolves from the equi- 
librium state. Substituting the equilibrium distribution 



function into Eq. (78) we get the stress tensor functional 



JV(£>*) 



5nv 



(81) 



which is proportional to the initial equilibrium pressure, 
-Po(0- The last step is to substitute the nonadiabatic 
"equation of state" , Eq. (81), into the force balance equa- 
tion of Eq. (76). This results in the following "hydrody- 
namic" equation of motion 



d 



U-m 



+ 



dp 



1 p dg aa 

2 dp 



= 0, 



(82) 



where we used the definition of the covariant divergence, 
Eq. (57), to compute the stress force, yfgP^ :vl in Eq. (76). 
We would like to outline that n (£) and Po(£) m Eq. (82) 
are the time independent initial density and pressure re- 
spectively. Equations (82) and (18) constitute a closed 
set of continuum mechanics equations which describe a 
long wavelength dynamics of a Fermi gas in the time- 
dependent Hartree approximation. Since the stress force 
in Eq. (82) depends only on the deformation tensor, g^, 
it is natural to interpret Eqs. (82) and (18) as a nonlin- 
ear elasticity theory of a Fermi gas. In the case of small 
deformations this theory reduces to the standard linear 
elasticity theory with a nonzero shear modulus. Indeed, 
in the linear regime Eq. (18) takes the form 



MM) 
dt 



(83) 



where u = x— £ is the displacement vector. The deforma- 
tion tensor reduces to the common linearized expression 



du^ du v 
~dp~~dp 



(84) 



Assuming for simplicity that the unperturbed state is 
homogeneous, and substituting Eqs. (83) and (84) into 
Eq. (82), we get the following equation of motion for the 
displacement vector 



d 2 u ll da^ dSU 
mn -—f- + no^rr- = 0. 



dt 2 dp 



dp 



(85) 



The linearized stress tensor <7 M „ takes the standard elastic 
form 



x v dUa 1 ( du v 



du v 2du a . 



where K = |Po an d £t = Po are the bulk modulus and 
the shear modulus of a Fermi gas respectively 19,20 ' 28 . 

Full nonlinear set of Eqs. (82), (18) is equivalent to 
the generalized collisionlcss hydrodynamics derived in 
Refs. 19,20 (see also Ref. 29). In fact, Eqs. (82), (18) 
and the generalized hydrodynamics of Refs. 19,20 corre- 
spond to the same theory in the Lagrangian and Eulerian 
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formulations respectively. An advantage of the present 
Lagrangian formulation is the explicit form of the stress 
tensor P^, Eq. (81). The Lagrangian point of view also 
gives a very clear microscopic picture of a fast collision- 
less dynamics. This is a kind of evolution of a many- 
body system with almost time-independent distribution 
of particles inside every moving and deforming fluid el- 
ement. Using the nonadiabatic equation of state in the 
Lagrangian frame, Eq. (81), we can easily recover the 
corresponding expression for the stress tensor P fiu (x,t) 
in the laboratory frame 



d£ a d£P ~ 

= .9^(x,t)v / 5(x,t) J Po(C(x,t)) 
Here ^^(x, t) is Cauchy's deformation tensor 1 



dxv dx v ' 



9 = 1/9 



(87) 



(88) 



One can check that P l _ lv (x, t) of Eq. (87) is a solution to 
the equation for the stress tensor derived in Refs. 19,20. 
The stress tensor P At!y (x, t), Eq. (87), enters the force bal- 
ance equation in the laboratory frame, Eq. (16). This 
is a hydrodynamic equation of motion in Eulerian de- 
scription. It is quite natural that P^ (x, t) depends on 
g M „(x, i) since Cauchy's tensor is a common characteris- 
tics of deformations in Eulerian picture. 

In contrast to the stress tensor in the Lagrangian 
frame, Eq. (81), the stress tensor of Eq. (87) is a highly 
nonlocal function. It is proportional to the pressure Po at 
the initial position of a fluid element which is currently at 
x (x is an independent variable in Eq. (16)). The local- 
ity of the stress force in Eq. (82) is a key property of the 
present Lagrangian formulation of a nonadiabatic contin- 
uum mechanics of a Fermi gas. This formulations should 
be much more convenient for applications to particular 
nonlinear problems. 



V. CONCLUSION 

We applied the idea of the Lagrangian description in 
continuum mechanics to the theory of nonequilibrium 
quantum many-body systems. Reformulation of the mi- 
croscopic many-body theory in terms of Lagrangian co- 
ordinates corresponds to the transformation to the local 
noninertial reference frame moving with the flow (the co- 
moving Lagrangian frame). This transformation allows 
to separate the convective motion of particles, which is 
a direct generalization of the common separation of the 
center-of-mass motion in homogeneous systems. The mo- 
tion of particles in the Lagrangian frame is influenced 
by the external forces and by generalized inertia forces. 
We have shown that the inertia forces can be described 
in purely geometric terms of Green's deformation tensor 
g^v and the skew-symmetric vorticity tensor F^. Ten- 
sors g^ u and enter equations of motion as an effec- 



tive metric tensor and an effective magnetic field respec- 
tively. Our results demonstrate a close relation of the 
many-body dynamics in Lagrangian frame to the quan- 
tum dynamics on curved manifolds. 

We also derived local conservation laws for the num- 
ber of particles and for momentum in the Lagrangian 
frame, and presented closed microscopic expressions for 
the stress tensor and for the corresponding stress force. 
The local momentum conservation law in the Lagrangian 
frame reduces to a zero force condition. The inertia forces 
exactly compensate the external force and the stress force 
in every point of the Lagrangian £-space. The net force, 
exerted on every fluid element, is exactly zero, which 
results in zero current density and a time-independent 
density distribution. This property is the main advan- 
tage of the Lagrangian description. It suggests one of 
the most promising application of our formalism, which 
is a new reformulation of TDDFT in a form similar to 
the static theory. Indeed the main practical problem of 
TDDFT is an inevitable strong nonlocality of exchange 
correlation potentials 4,7,8 . The physical reason for this 
is just the nonadiabatic motion of fluid elements. When 
time is flowing, new and new fluid elements arrive at 
a given point x, and bring an information about sur- 
rounding space, producing the above nonlocality. Using 
our reformulation of the many-body theory as a basis for 
TDDFT one can completely remove the very source on 
the nonlocality, which is of extreme practical importance. 
In this paper we did not touch these questions since it 
required an extended special consideration. A detailed 
formulation of TDDFT in the Lagrangian frame will be 
presented in the next paper of this series 30 . 

In this paper we also considered two illustrative exam- 
ples of application. The most interesting of them is the 
description a nonlinear semiclassical dynamics of a colli- 
sionless Fermi gas. We have shown that the full problem 
can be separated into two independent parts. The first 
one is the solution of a universal kinetic problem, which 
defines the stress tensor as a universal functional of g^v 
and Ff,, v . This stress tensor is used as an input for the 
second "hydrodynamic" part of the problem, determin- 
ing the dynamics of the velocity vector. This separa- 
tion of the initial many-body problem can be viewed as 
a particular realization of TDDFT in the hydrodynamic 
formulation 4 ' 22 . In the case of a fast long wavelength 
dynamics (similar to that for plasma oscillations), the 
universal kinetic problem can be solved explicitly. The 
solution is extremely simple - the Wigner function in the 
Lagrangian frame is time- independent. The correspond- 
ing "hydrodynamic" problem also can be formulated in 
the explicit form. It reduces to a closed nonlinear elas- 
ticity theory of a Fermi gas. This elasticity theory is, 
in fact, the Lagrangian formulation of the generalized 
hydrodynamics derived in Refs. 19,20. The generalized 
hydrodynamics proved to be useful in the description of a 
small-amplitude collective dynamics of an electron gas. It 
gives the correct (consistent with the kinetic treatment) 
dispersion of plasma waves in homogeneous systems 19 ' 20 
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and recovers the exact dispersion of the edge modes in a 
confined geometry 31 . The results for the standing plasma 
waves in a parabolically trapped electron gas are also 
quite reasonable 32,33 . The Lagrangian "elasticity theory" 
of a Fermi gas, derived in Sec. IVB, is structurally much 
more simple than the Eulerian formulation of Refs. 19,20. 
Therefore, we believe that it should provide a good ba- 
sis for the description of nonlinear dynamical effects in 
inhomogeneous many-electron systems. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft under Grant No. PA 516/2-3. 



APPENDIX A: THE DIVERGENCE 
REPRESENTATION OF THE INTERACTION 
STRESS FORCE 

By definition the infraction stress force F lnt is the com- 
mutator of the current operator j and the interaction 
Hamiltonian W , 



F int (x)=zm[j(x),IF]. 



(Al) 



In the coordinate representation operators j and W are 
defined as follows 

JW = iE[p^( x - x -)]+- ( A2 ) 

i 

W=^w(x i ,x j ), (A3) 

where i and j label particles. Calculating the commuta- 
tor of Eq. (Al) we get the force in the following form 

^^(Xj.Xj) 



1,3 



+S(x-Xj)- 



(A4) 



If the interaction potential satisfies the Newton's third 
law, 



dw(xi,Xj) dw(xi,Xj) 



(A5) 



Eq. (A4) takes the form 

F int (x) = \ £ [5(x - Xi) - 6(x - x,)] dw ^\ 



(A6) 

The difference of delta-functions in Eq. (A6) can be trans- 
formed as follows 



6(x — Xi) — S(x — Xj) 



1 — e ( x i _x i)7^ 



S(x - Xi) 



1 

= -(x i -x J )^ J d\e^-^6(x-x<) 
o 

i 

d f 

= -^( x » - Xj) J dXS(x - x, - X(xj - x^). (A7) 

o 

Inserting Eq. (A7) into Eq. (A6) we get the final repre- 
sentation for the interaction stress force 



d 



(A8) 



where W M „(x) is the interaction stress tensor operator 
l 

W^(x) = ~^2 J d\ S(x-Xi - X(xj -xt)) 



x«-<)- 



The A-integration in Eq. (A9) is along the line that con- 
nects two interacting particles. 
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